Two-loop renormalization of the effective field theory 

of a static quark 



D. J. BROADHURST, A. G. GROZIN* 



Physics Department, Open University, Milton Keynes MK7 6AA, UK 



Abstract We give a recurrence relation for two-loop integrals encountered in the 
effective held theory of an inhnitely heavy quark, Q, interacting with gluons and A^l 
massless quarks, g, from which we obtain exact two-loop results, in any dimension 
and covariant gauge, for the propagator of Q and the vertex function of the heavy- 
light current J = QTq, at zero g-momentum. The anomalous dimension of the Q- 
held agrees with the recent result of Broadhurst, Gray and Schilcher. The anomalous 
dimension of the current is 

^ dlogZj a, ( , /127 + 56C(2)-10iVL\ as , 2 

7J = -n = 1 1 + ^7; ^ 0{as 

dlog^ vr I V 72 J TT ^ 

which gives the new two-loop correction to the result of Voloshin and Shifman. 
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1 Introduction 



Recently an interesting new approach to QCD problems involving a heavy quark 
was proposed, namely the effective field theory (EFT) of a static quark [1-5]. (See 
also the review [^].) To leading order in the 1/m expansion, the EFT lagrangian 

L = QliDoQo + %iPqo - \Gl^^Gl^^ - (d^A^^^) ^/2ao + (ghost term) (1) 

where qq and Aq^ are the bare light-quark and gluon fields of conventional QCD and 
Qo is the bare static-quark field, which is a two-component spinor. (When used as 
a 4-component spinor, its lower components are assumed to vanish.) The free-field 
propagator of Qo is l/<^ (equivalently, (1 + ^q)/2uj, in 4-component form) where 
u) = pq — ra \s the energy by which it is off-shell. Its coupling to gluons is given by 
the vertex igQV^f^^ where = (1,0) is its velocity. The lagrangian (||) possess an 
SU(2) static-quark spin symmetry Q. 

In this paper we consider the two- loop renormalization of EFT, using dimen- 
sional regularization in D = 4 — 2e dimensions. In the MS scheme the bare quantities 
in (|l|) are related to the corresponding renormalized quantities by 

Qo = n^'zl'^Q, qo = -^-'Zl'\, Ag^ = U'^ZI'^AI, go = U'zH^g, ao = Z^a 

where Z2, Z^ and Za are the same as in QCD [^, with A'^l light-quark flavours, 
and fi^ = /U^e'^/47r. Here we calculate, to two loops, the static-quark wave-function 
renormalization constant Z2 and the renormalization constant Zj of the heavy-light 
current Jo = Qo^^^o = ZjJ, where P is some (irrelevant) gamma matrix. For 
Z2 we obtain the same result as Broadhurst, Gray and Schilcher who have 
independently extracted it from the gauge-invariant singularities of on-shell wave- 
function renormalization of a finite-mass quark P]. ^From Zj we obtain the two- 



loop correction to the leading-order result [10, 11| for the anomalous dimension 
7j = d log Zj /d log fj,. This is needed for the renormalization-group analysis of the 
extrapolations of lattice computations of heavy-meson decay constants jl^]. 

Renormalization of EFT corresponds closely to renormalization of the Wilson 
line [13-17], which coincides with the static-quark propagator in coordinate space. 
In particular the lagrangian (1) is essentially the same as that proposed in [U], for 
investigating Wilson lines, and used in the two-loop calculations of |jl^, ^]. In 
Aoyama calculated Z2 in the Feynman gauge, at Ni^ = 0; we disagree with his result. 
In [17] Korchemsky and Radyushkin calculated the two- loop anomalous dimension 



of a heavy-heavy current, which depends on the relative velocity of initial and final 
quarks; the corresponding one- loop result [Q has recently been rediscovered in 
EFT M. 



2 One- and two-loop integrals 

The static-quark bare self-energy term —[T,{uj) is given, to two loops, by the diagrams 
of Fig 1, where the blob in Fig la includes the light-quark, gluon and ghost loops in 
the gluon propagator. We take the light quarks as massless and easily evaluate the 
master one-loop EFT integral 
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using the parameterization 

1 r(a+p) y^'^dy 



A'^PP T{a)T{p) Jo {A + 

where we adopt the convention that massless hnes have greek indices and static Unes 
have roman ones. We find that 

Ti2a+p-D)T{D/2-a) 

'^"'^^ = rW¥) 

which enables us to evaluate diagrams of the form of Fig la. 

The master two-loop EFT integral, for diagrams of the form of Fig lb, is 



l\a/_l\/3/ _1 



LJ \^ IjJ 



k"^ J \P J \{k-lYJ \uj + koJ \uj + k 
7r^(-2..)2(^-"-^-^)/(a,/3,7,p,g) 



In degenerate cases one uses (Q) and the corresponding result |15] for the master 
one-loop diagram of massless QCD 

r{a + P- D/2)T{D/2 - a)T{D/2 - /?) 



T{a)T{f3)T{D - a - (5) 



to obtain 



I(a,/3,0,p,g) = I{a,p)I{(3,q) (3) 
/(a,0,7,p,g) = I{a,2j+p + q- D)I{j,q) (4) 
/(a,/3,7,P,0) = I(a + /3 + 7-^/2,p)G(/3,7) (5) 

which suffice for diagrams of the form of Fig Ic and Fig Id, since the integrand of 
Fig Ic has only 4 factors in the denominator and the integrand of Fig Id can be 
expressed as the sum of three such terms after multiplying by 

1 = -{{UJ + ko) + {UJ + lo) -{LJ + ko + lo)} 
to 

This is an enormous simplification in comparison with on-shell massive QCD calcu- 
lations [^, ^, where this last diagram was the most difficult. 

In general we use a recurrence relation to reduce the non-degenerate case of the 
master integral to a sum of known degenerate cases. Using the method of integration 
by parts [^, |l8|, we have found that 

{D-a-2-f-p-q+l)I = {QA+(r- - B") + (2(D - q - /3 - 7) - p - g + 1)Q-} / 

(6) 

where, for example. A"*" raises a and Q lowers q. Again this is simpler than the 
triangle relation ||l^ of massless QCD, because no more than three terms appear 
and no static-line index is increased. 

Now it is a matter of programming to reduce every one-loop integral to a rational 
function of D = 4 — 2e times 

Fi = (-2cj)-2T(-e)F(l + 2e)/(4^)^/2 (7) 
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and every two-loop integral from Fig 1 to a linear combination of Vf and 



r2 = (-2a;)-^T2(-e)r(l + 4e)/(4^)^ (8) 
= (l + 4C(2)e2-16C(3)e3)r2 + 0(e2) (9) 

with coefficients which are rational functions oi D. We used REDUCE |T|] on a 
Vax and checked results with Mathematica [^] on a Macintosh. The programming 
effort and run times were far less than in the related QCD calculations ^. 



3 Anomalous dimension of the static-quark field 

Using these methods, it is straightforward to compute the two- loop corrections to 
the bare propagator Sq{uj). We find that 

ujSoiu;) = ^— - = 1 + ZiEiCf^o + (^aiCaE? 

+ (^A2Ca + Zf2Cf + ZL2rFiVL)r2} Cp^O + ^(^o) (10) 



with coefficients 
Zi = ao - 



D 



Za 



D - 
Zi 



{D - 2f{D - 5) (D^ -AD + 5)Zi (D^ - 9D + 16) Zf 

~ 2(2D-7)(L>-3)3(L»-6) ^ 2(L» - 3)2(1) - 6) ~ 8{D - 3){D - 6) 

2{D - 2) 



Zf2 
Zh2 



{2D -7){D -3){D -Q) 



and colour factors Ca = N, Cp = {N"^ — 1)/2N, Tp = ^, for an SU(A^) gauge group, 
or Ca = 0, Cp = Tp = 1 for a U(l) gauge group. 

__ ~i/2 

We now renormalize the bare static-quark field Qq = ji ^Z^ Q using the mini- 
mal wave-function renormalization constant Z2 which makes the renormalized prop- 
agator S{u) = S{){uj)/Z2 finite. (Tildes are used here to distinguish EFT renormal- 
ization constants from those of QCD.) We work in the MS scheme, in which the bare 
coupling and gauge parameter are renormalized by 0] 

^° - Z^ = l-^(^4cA-\TFN^+0{al) (11) 



g'^p?' 47re V 3 3 

5! = Z3.1 + i^(H^C^-lrp%)+0K=) (12) 

whereas = g'^/A-K. We posit a minimal form Z2 = l + (Cii(a)/e)as/7r + (C22(a)/e^ + 
C2i{a) / e)a1 / tt"^ + 0{a^) and determine its coefficients from (p^. The renormaliz- 
ability of the theory requires that this constant make Sq{ijj)/ Z2 finite for all a;, which 
provides a highly non-trivial check on our result. Then, using the jj. dependence of as 



and o, implied by (|ll| , [12D , we obtain the anomalous dimension 7p = d log Z2 /d log 



of the static Q field. As a check on the whole procedure we mirror it, step by step. 
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by that for a massless quark, to obtain the corresponding QCD quantity 7f for the 
massless q field from the recurrence relations of Our result is 

' 27r I 48 16 12 J 7r2 ^ ^ ^ 

where the latter form was independently obtained in |^] from the gauge-invariant 
singularities of on-shell wave-function renormalization of a finite-mass quark. Note 
that (^) vanishes in QED in the renormalized Yennie gauge [^] o = 3, for which 



there is no 'infrared catastrophe' [22|. From ( p~3| ) we find that 7f agrees with the 
QCD result of Egorian and Tarasov ||2^ (misquoted in which contains sign 
errors). 

Setting a = 1 and A^L = in (||), we discover an error in [|16[, where the 



renormalization of the Wilson line was calculated (neglecting light-quark loops) in 
the Feynman gauge (where Fig lb gives no contribution). The Oi^/e^ terms in |[l^ 
are inconsistent with the renormalization group and the ag/e terms include the 
contributions of Figs lc,d with the wrong sign. This incorrect result was used 
in [p^ j, where the authors also omitted a factor of 2 from the two- loop term in 
the anomalous dimension. In view of these approximations and errors the two-loop 



analysis of |24| needs reconsideration. 



4 Anomalous dimension of the heavy-light current 

The simplest method for finding the anomalous dimension 7j of the current Jq = 
QqTqq is to consider the bare vertex function with zero light-quark momentum. The 
diagrams of Fig 2 then differ from those of Fig 1 only in their numerator structure 
and the indices of their denominators. The effect is to multiply the gamma matrix 
r by a factor Tq{uj), which is gauge-dependent and more complicated than (p^). 
But it becomes much simpler when one multiplies by the external static-quark bare 
propagator So{ui). We find that 

Vo{lo) = LoSo{io)ro{ij) = i + ViriCFgl + {{VAiCA + VFiCF)rl 

+ {Va2Ca + Vf2Cf + VL2TFNi^)r2} Cpg^ + 0{gl) (14) 

is gauge-invariant. The coefficients are 

D - 1 



Vi 
Vai 
Vfi 

Va2 
Vf2 
Vl2 



D-3 

-4D-1 
2{D - 3)2(D -4) 

D-2 
{D-3){D-4) 

2D^ - 22D^ + 891)2 _ ^^^j^ _^ 
~ 2(2i:> - 7){D - 3)2(1) - A){D - 6) 

(Z)2 -bD + 2)D 
~2(L»-3)2(Z)-4)(L>-6) 
2(D - 2) 



{2D - 7){D - 3){D - 6) 
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independently of the spin structure of the current J. 

The renormahzation constant Zj = Z^Z^"^ Z^"^ may be calculated from the re- 
quirement that the renormalized vertex function r(a;) = Tq{ijo)/Z-p be finite. Equiv- 
alently, but more conveniently, we have Zj = ZyZ2 Z2 , where Z\ = Zr^2 
makes V{oj) = Vq{lo)/Z\i finite. This gives two highly non-trivial checks on our re- 
sults. First, the ot1/e^ terms in (|l^ must be such that division by a constant makes 
V{uj) finite for all oj. Second, the anomalous dimension 7v = d log /d log /i must 
be gauge-invariant, to ensure that 7j = 7v — ^7f + ^7f is also gauge-invariant. As 
a third check, we mirrored this EFT calculation, step by step, by the QCD calcula- 
tion of 7gq, obtaining agreement with Tarrach |2^ ]. These checks give us considerable 
confidence in our final result 

= _^^_/1^Ca--Cf- ArFiVL-aCA-CF)C(2)l^ + 0(a3) 

= l^, + {iCA + iCF + (iCA-CF)c(2)}^ + 0(a.^) (15) 

whose second form we have also verified from the quark-condensate contributions to 
the correlator of J |26]. 



5 Summary and conclusions 

We have given a method for calculating one-scale one- and two-loop EFT integrals. 
The master two-loop integral can be reduced by the recurrence relation (P) to the 
degenerate cases (|^,^J5|) , which yield the structures . We have used this method 
to calculate the static-quark propagator and a particular case of the heavy-light 
vertex function, obtaining (^Jl^) in an arbitrary spacetime dimension and in any 
covariant gauge. From these follow the anomalous dimensions of the static-quark 
propagator (13) and the heavy-light current ([T5[). The first agrees with Q and 
contradicts an old result for the Wilson line ||l^ . The second is also confirmed in 1 2C ] . 



It gives the radiative correction to the one-loop result of |10, 11] and is necessary for 
matching QCD and EFT currents, with account of the radiative corrections of [|^. 
The lack of it in [13] limited the accuracy of extrapolation of lattice computation. 



Novel features of EFT calculations include: the trivial simplification of diagrams 
such as Figs Id and 2i, which are the most difficult in on-shell massive QCD calcu- 
lations |§, an algorithm more efficient than the triangle relation |l^] of massless 
QCD; the appearance of ("(2) = 7r^/6 in counterterms, via the relation ^ between 
the basic integral structures, in contrast to the situation in QCD, where C(3) is the 
first term to appear in the expansion. 
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